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_ Abstract 

^ I A results of numerical procedure for construction of confidence intervals for parameter of 

'^ I Poisson distribution for signal in the presence of background which has Poisson 
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distribution with known value of parameter are presented. It is shown that the described 
procedure has both Bayesian and frequentist interpretations. 
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g : I. INTRODUCTION 

o 

O ' In paper illl the unified approach to the construction of confidence intervals and confi- 

>~>' dence limits for a signal with a background presence, in particular for Poisson distributions, 

R-'l has been proposed. The method is widely used for the presentation of physical results 

k>( ; though a number of investigators criticize this approach [^ 

C^ I In present paper we use a simple method for construction of confidence intervals for 

parameter of Poisson distribution for signal in the presence of background which has Poisson 
distribution with known value of parameter. This method is based on the statement M that 
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the true value of parameter of the Poisson distribution in the case of observed number of 
events x has a Gamma distribution. In contrast to the approach proposed in [^, the width 
of confidence intervals in the case of a; = is independent on the value of the parameter of 
the background distribution. The described procedure has both Bayesian and frequentist 
interpretations. 

In Section 2 the method of construction of confidence intervals for parameter of Poisson 
distribution for signal in the presence of background which has Poisson distribution with 
known value of parameter is described. The results of confidence intervals construction and 
their comparison with the results of unified approach are also given in the Section 2. The 
main results of this note are formulated in the Conclusion. 

II. THE METHOD OF CONSTRUCTION OF CONFIDENCE INTERVALS 

Assume that in the experiment with the fixed integral luminosity (i.e. a process under 
study may be considered as a homogeneous process during given time) the x events of some 
Poisson process were observed. It means that we have an experimental estimation \{x) 
of the parameter A of Poisson distribution. We have to construct a confidence interval 
(Ai(x), A2(a;)), covering the true value of the parameter A of the distribution under study 
with confidence level 1 — a, where a is a significance level. It is known from the theory of 
statistics 0, that the mean value of a sample of data is an unbiased estimation of mean 
of distribution under study. In our case the sample consists of one observation x. For the 
discrete Poisson distribution the mean coincides with the estimation of parameter value, 
i.e. A = a: in our case. As it is shown in ref [Q the true value of parameter A has Gamma 
distribution Pi £+i, where the scale parameter is equal to 1 and the shape parameter is equal 
to X + 1 (see Fig.l), i.e. 

PiX\x) = P{x\X) = ^e-''. (2.1) 

x\ 

Note that formula (2.1) results from the Bayesian formula H 



P{X\x)P{x) = P{x\X)P{X) (2.2) 

in the assumption that all possible values of parameter A have equal probability, i.e. P{X) = 
const. In this assumption the probability that unknown parameter A obeys the inequalities 
Ai < A < A2 is given by evident Bayesian formula 

P(Ai < A < Asix) = P(Ai < A|x) - P(A2 < A|x) = / ' P{X\x)dX, (2.3) 

P{X, < X\x) = J^ P{X\x)dX, 

where P{X\x) is determined by formula (2.1). 

Formula (2.3) has also well defined frequentist meaning. Using the identity 

00 \ip-Xl i-\2 \^p-^ ^ AV~'^2 



i=x+l ^- "^^1 ^- i=0 ^• 

one can rewrite formula (2.3) as 

P(Ai < A < X2\x) = l-P{n<x\X2)-P{n> x\Xi)=P{n<x\Xi)-P{n<x\X2), 

(2.5) 

X \n —\ 00 \n —\ 

where P{n < x\X) = ^ j — and P{n > x\X) = ^ j — . 

n=0 ''^- n=x+l ""• 

The right hand side of formula (2.5) has well defined frequentist meaning and in fact it 
is one of the possible definitions of the confidence interval in frequentist approach Q As an 
example, such type the shortest 90% CL confidence interval in case of observed number of 
events x = 4 is shown in Fig. 2. 

For instance, for the case A2 = cxd formula (2.5) takes the form 

/ P{X\x)dX = P{Xi<X\x) = Y.P{n\Xi) = P{n<x\Xi), (2.6) 

"'^1 n=0 

which has evident frequentist meaning too. 



^See, however, ref. 



Let us consider the Poisson distribution with two components: signal component with 
a parameter A^ and background component with a parameter A;,, where A;, is known. To 
construct confidence intervals for parameter A^ of signal in the case of observed value x we 
must find the distribution P{Xs\x). 

At first let us consider the simplest case x = s + b = 1. Here s is a number of signal 
events and 6 is a number of background events among observed x events. 

The b can be equal to and to 1. We know that the b is equal to with probability 

Po = P{b = 0) = ^e-^'^ = e-^" (2.7) 

and the b is equal to 1 with probability 

Pi = P{b = 1) = ^^e-"" 

Correspondingly, P(b = 0\x = 1) = P(s = llx = 1 

Po + Pi 

P(s = 0|x = 1) = ^^ . 
Po+Pi 

It means that distribution of P{Xs\x = 1) is equal to sum of distributions 

p{s = i\x = i)ri,2 + P{s = o\x = i)ri,i = ^^ri,2 + ^^ri,i, (2.9) 

Po + Pi Po + Pi 

where Fi i is Gamma distribution with probability density P(As|s = 0) = e"'^" and ri2 is 
Gamma distribution with probability density P(As|s = 1) = Ase~'^^ As a result we have 

P(A,|x=l) = ^i±^e-^^ (2.10) 

i + Af, 

Using formula (2.10) for P{Xs\x = 1) and formula (2.5) we construct the shortest confi- 
dence interval of any confidence level in a trivial way. 

In this manner we can construct the distribution of P{Xs\x) for any values of x and A;,. 
As a result we have obtained the known formula fj. 



A5e-^^ 




(2.8) 


^1) = : 


Po 


- and P{b = 1 X = 1) = 



•fThe formula (2.11) has been derived earlier in ref. H (formula 5.88) in the framework of Bayesian 
approach and in ref. |^]. We thank Prof. D'Agostini for correspondence. 



P{K\x) = ^^^4^e-^ (2.11) 
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a:!E- 



i=0 

The numerical results for the confidence intervals and for comparison the results of 
paper |jl| are presented in Table 1 and Table 2. 

It should be noted that in our approach the dependence of the width of confidence 
intervals for parameter A^ on the value of A^ in the case £ = is absent. For x = the 
method proposed in ref. |T^ also gives a 90% upper limit independent of A^. This dependence 



is absent also in Bayesian approach 0,11 



III. CONCLUSION 

The results of construction of frequentist confidence intervals for the parameter A^ of 
Poisson distribution for the signal in the presence of background with known value of pa- 
rameter Afe are presented. It is shown that the described procedure has both Bayesian and 
frequentist interpretations. 
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FIGURES 




FIG. 1. The behaviour of the probabihty density of true value of parameter A for Poisson 



A^ 



-A • 



distribution in case of x observed events versus A and x. Here fix, A) = — -e is both a Poisson 

x\ 

distribution with parameter A along the axis x and a Gamma distribution with a shape parameter 
X + 1 and a scale parameter 1 along the axis A. 




FIG. 2. The Poisson distributions f{x, A) for A's determined by the confidence hmits Ai = 1.51 
and A2 = 8.36 in case of observed number of events x = 4 are shown. The probabiUty density of 
Gamma distribution with scale parameter a = 1 and shape parameter x = x = 4 is shown inside 
this confidence interval. 



TABLES 
TABLE L 90% C.L. intervals for the Poisson signal mean As, for total events observed x, for 

known mean background A^ ranging from to 4. A comparison between results of ref.[l] and results 

from present note. 



x\X, 


0.0 rcf.[l] 


0.0 


1.0 rcf.[l] 


1.0 


2.0 rcf.[l] 


2.0 


3.0rcf.[l] 


3.0 


4.0ref.[l] 


4.0 





0.00, 2.44 


0.00, 2.30 


0.00, 1.61 


0.00, 2.30 


0.00, 1.26 


0.00, 2.30 


0.00, 1.08 


0.00, 2.30 


0.00, 


1.01 


0.00, 2.30 


1 


0.11, 4.36 


0.09, 3.93 


0.00, 3.36 


0.00, 3.27 


0.00, 2.53 


0.00, 3.00 


0.00, 1.88 


0.00, 2.84 


0.00, 


1.39 


0.00, 2.74 


2 


0.53, 5.91 


0.44, 5.48 


0.00, 4.91 


0.00, 4.44 


0.00, 3.91 


0.00, 3.88 


0.00, 3.04 


0.00, 3.53 


0.00, 


2.33 


0.00, 3.29 


3 


1.10, 7.42 


0.93, 6.94 


0.10, 6.42 


0.00, 5.71 


0.00, 5.42 


0.00, 4.93 


0.00, 4.42 


0.00, 4.36 


0.00, 


3.53 


0.00, 3.97 


4 


1.47, 8.60 


1.51, 8.36 


0.74, 7.60 


0.51, 7.29 


0.00, 6.60 


0.00, 6.09 


0.00, 5.60 


0.00, 5.35 


0.00, 


4.60 


0.00, 4.78 


5 


1.84, 9.99 


2.12, 9.71 


1.25, 8.99 


1.15, 8.73 


0.43, 7.99 


0.20, 7.47 


0.00, 6.99 


0.00, 6.44 


0.00, 


5.99 


0.00, 5.72 


6 


2.21,11.47 


2.78,11.05 


1.61,10.47 


1.79,10.07 


1.08, 9.47 


0.83, 9.01 


0.15, 8.47 


0.00, 7.60 


0.00, 


7.47 


0.00, 6.76 


7 


3.56,12.53 


3.47,12.38 


2.56,11.53 


2.47,11.38 


1.59,10.53 


1.49,10.37 


0.89, 9.53 


0.57, 9.20 


0.00, 


8.53 


0.00, 7.88 


8 


3.96,13.99 


4.16,13.65 


2.96,12.99 


3.18,12.68 


2.14,11.99 


2.20,11.69 


1.51,10.99 


1.21,10.60 


0.66, 


9.99 


0.34, 9.33 


9 


4.36,15.30 


4.91,14.95 


3.36,14.30 


3.91,13.96 


2.53,13.30 


2.90,12.94 


1.88,12.30 


1.92,11.94 


1.33,11.30 


0.97,10.81 


10 


5.50,16.50 


5.64,16.21 


4.50,15.50 


4.66,15.22 


3.50,14.50 


3.66,14.22 


2.63,13.50 


2.64,13.21 


1.94,12.50 


1.67,12.16 


20 


13.55,28.52 


13.50,28.33 


12.55,27.52 


12.53,27.34 


11.55,26.52 


11.53,26.34 


10.55,25.52 


10.53,25.34 


9.55,24.52 


9.53,24.34 



TABLE IL 90% C.L. intervals for the Poisson signal mean A^, for total events observed x, for 
known mean background Aj, ranging from 6 to 15. A comparison between results of ref.[l] and 
results from present note. 



x\Xt 


6.0rcf.[l] 


6.0 


8.0rcf.[l] 


8.0 


10.0 rcf.[l] 


10.0 


12.0 rcf.[l] 


12.0 


15.0rcf.[l] 


15.0 





0.00, 


0.97 


0.00, 


2.30 


0.00, 


0.94 


0.00, 


2.30 


0.00, 


0.93 


0.00, 


2.30 


0.00, 


0.92 


0.00, 


2.30 


0.00, 


0.92 


0.00, 2.30 


1 


0.00, 


1.14 


0.00, 


2.63 


0.00, 


1.07 


0.00, 


2.56 


0.00, 


1.03 


0.00, 


2.51 


0.00, 


1.00 


0.00, 


2.48 


0.00, 


0.98 


0.00, 2.45 


2 


0.00, 


1.57 


0.00, 


3.01 


0.00, 


1.27 


0.00, 


2.85 


0.00, 


1.15 


0.00, 


2.75 


0.00, 


1.09 


0.00, 


2.68 


0.00, 


1.05 


0.00, 2.61 


3 


0.00, 


2.14 


0.00, 


3.48 


0.00, 


1.49 


0.00, 


3.20 


0.00, 


1.29 


0.00, 


3.02 


0.00, 


1.21 


0.00, 


2.91 


0.00, 


1.14 


0.00, 2.78 


4 


0.00, 


2.83 


0.00, 


4.04 


0.00, 


1.98 


0.00, 


3.61 


0.00, 


1.57 


0.00, 


3.34 


0.00, 


1.37 


0.00, 


3.16 


0.00, 


1.24 


0.00, 2.98 


5 


0.00, 


4.07 


0.00, 


4.71 


0.00, 


2.60 


0.00, 


4.10 


0.00, 


1.85 


0.00, 


3.72 


0.00, 


1.58 


0.00, 


3.46 


0.00, 


1.32 


0.00, 3.20 


6 


0.00, 


5.47 


0.00, 


5.49 


0.00, 


3.73 


0.00, 


4.67 


0.00, 


2.40 


0.00, 


4.15 


0.00, 


1.86 


0.00, 


3.80 


0.00, 


1.47 


0.00, 3.46 


7 


0.00, 


6.53 


0.00, 


6.38 


0.00, 


4.58 


0.00, 


5.34 


0.00, 


3.26 


0.00, 


4.65 


0.00, 


2.23 


0.00, 


4.19 


0.00, 


1.69 


0.00, 3.74 


8 


0.00, 


7.99 


0.00, 


7.35 


0.00, 


5.99 


0.00, 


6.10 


0.00, 


4.22 


0.00, 


5.23 


0.00, 


2.83 


0.00, 


4.64 


0.00, 


1.95 


0.00, 4.06 


9 


0.00, 


9.30 


0.00, 


8.41 


0.00, 


7.30 


0.00, 


6.95 


0.00, 


5.30 


0.00, 


5.89 


0.00, 


3.93 


0.00, 


5.15 


0.00, 


2.45 


0.00, 4.42 


10 


0.22,10.50 


0.02, 


9.53 


0.00, 


8.50 


0.00, 


7.88 


0.00, 


6.50 


0.00, 


6.63 


0.00, 


4.71 


0.00, 


5.73 


0.00, 


3.00 


0.00, 4.83 


20 


7.55,22.52 


7.53,22.34 


5.55,20.52 


5.53,20.34 


3.55,18.52 


3.55,18.30 


2.23,16.52 


1.70,16.08 


0.00,13.52 


0.00,12.31 



